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We argue that dark matter particles which have strong interactions with the Standard Model 
particles are not excluded by current astrophysical constraints. These dark matter particles have 
unique signatures at colliders; instead of missing energy, the dark matter particles produce jets. We 
propose a new search strategy for such strongly interacting particles by looking for a signal of two 
trackless jets. We show that suitable cuts can plausibly allow us to find these signals at the LHC 
even in early data. 
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Introduction. Many different astrophysical and cos- 
mological measurements have definitively established the 
existence and abundance of dark matter (DM) [l[ . How- 
ever, these observations do little to establish the masses 
and interaction strengths of these DM particles. A pri- 
ori, the DM particles only need to have gravitational in- 
teractions; however, much stronger interaction strengths 
are not excluded. Similarly, DM modelspropose masses 
ranging from sub-ev 0] to several TeV Q. 

Perhaps the most popular candidate for DM is a 
weakly interacting massive particle (WIMP), interacting 
via weak interactions with the Standard Model (SM) par- 
ticles [3| . Many models of physics beyond the SM (such 
as supersymmetry and technicolor) naturally have sta- 
ble WIMP candidates for DM. Furthermore, WIMPs of 
a mass close to a few hundred GeV naturally obtain the 
correct observed relic density through thermal interac- 
tions. Motivated by these considerations, many experi- 
ments have searched for WIMPs through direct detection 
of WIMP scattering off nuclei. However ,no unambiguous 
evidence of DM has been found so far [j, Q . 

WIMPs can also be searched for at colliders like the 
Tevatron and the Large Hadron Collider (LHC). Because 
of its weakly-interacting strength with ordinary matter, 
the WIMPs escape the detector and are undetected. This 
can lead to events where the visible particles have a 
nonzero total transverse momentum (balanced by the es- 
caping DM particles). Both the LHC and Tevatron have 
searched for such events with a large missing transverse 
momentum, and have seen no deviations from the predic- 
tions of the SM [6j, [Z| so far. This has placed significant 
constraints on many beyond-the-standard models with 
stable WIMP particles @-Q3. 

In view of these constraints, it may be worthwhile to 
reconsider alternative models of DM. It is particularly in- 
teresting to consider models where the DM is a strongly 
interacting massive particle (SIMP) [llj ]. Indeed, DM 
with strong interactions has features which may enable 
it to be undetected in searches for WIMPs. For exam- 
ple, SIMPs may be stopped in the atmosphere and earth, 
and cannot reach the detectors hundreds of meters under- 
ground. It is therefore worthwhile to see if some portion 
of SIMP parameter space is still viable. 



If the SIMPs have masses below a TeV, the Teva- 
tron and LHC can copiously produce them in pairs (if 
their cosmological stability is protected by a Z2 symme- 
try). Unlike WIMPs, SIMPs could behave like a neu- 
tron and deposit energy or even stop in the calorimeters. 
Therefore, for a wide range of parameter space, SIMPs 
produce two jets without missing energy, and the cur- 
rent DM searches with missing energy are insensitive to 
SIMPs. Compared to the QCD jets, the jets from SIMPs 
have zero tracks and less electromagnetic activity in the 
electromagnetic calorimeter and therefore can be distin- 
guished from the QCD dijet backgrounds. We will call 
these jets "dark matter jets". 

We will begin by considering various constraints which 
have already been placed on the SIMP parameter space. 
These come from many different varieties of astrophysical 
constraints, and we review these in the following section. 
We find that there are significant regions of parameter 
space which are still allowed. We then discuss the strate- 
gies and the prospects of probing these regions of param- 
eter space at the LHC. 

Constraints on SIMPs. For the purposes of dis- 
cussing constraints on SIMP models, it will often be con- 
venient to consider a model where the SM is extended by 
including a single DM particle x with a mass m x . An- 
other important parameter is the total cross section of 
DM scattering off nucleons a xp . Various constraints can 
be applied to the m x -a xp plane. We now consider the 
most important constraints in turn. 

Direct Detection Searches: Some of the most robust 
bounds on the SIMP parameter space come from di- 
rect detection searches. Ground-based experiments like 
CDMS and XENON place stringent bounds at some- 
what lower cross sections. At higher cross sections, as 
already mentioned, the DM particles are stopped by 
the atmosphere and the ground-based experiments be- 
come insensitive to them. However, high-altitude de- 
tectors (balloon-, satellite- or rocket-based experiments) 
can search for interactions above the atmosphere, and 
these eliminate large swaths of the parameter space. 
In particular, the RRS (balloon-based silicon semicon- 
ductor detector) [l2| and X-ray Quantum Calorimetry 
(XQC) 0, [3 experiments exclude SIMPs over a large 
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region of parameter space. We show the excluded regions 
at 90% C.L. in the enclosed regions (the blue and green 
lines) of Fig. [U 

Earth Heating Bounds: A particularly interesting 
bound comes from constraints on heat emitted from the 
Earth's core. SIMPs can be captured by the Earth and 
accumulate over time at the Earth's core. If SIMPs can 
self- annihilate into SM particles, the Earth's heat flow, 
which can be measured by placing detectors in deep un- 
derground shafts, would be modified. Mack and Bea- 
com [l5[ showed that this constraint effectively closes all 
of SIMP parameter space. However, this bound does not 
apply if SIMP self-annihilation is absent; for example 
asymmetric SIMP DM would evade this bound fl6| . 

Neutron Stars and Black Holes: Another interesting 
bound can be imposed for scalar DM; these can be col- 
lected in the cores of neutron stars and cause neutron 
stars to collapse to black holes. This excludes scalar 
SIMPs to very low masses [13, EH]- However, if the DM 
is fermionic, then this constraint is avoided. 

Cosmic Rays: Another important constraint comes 
from cosmic rays (l9| . Protons in cosmic rays can scat- 
ter off DM particles and produce neutral pions, which 
would decay to photons and be visible in gamma ray 
telescopes. Using this fact, the authors of [l9| placed a 
bound on DM-nuclcon interactions. However, these con- 
straints depend on many factors; most importantly they 
assume a form of the DM density near the galactic core. 
Since we are now dealing with a nonstandard form of 
DM with relatively strong interaction with baryons, these 
densities may be significantly modified, making the cos- 
mic ray constraints somewhat uncertain. In any case, we 
show the upper exclusion limit in the black line of Fig. [T] 
by requiring the cosmic ray contributions from DM to 
not exceed the observed flux. 

CMB and LSS: The large scattering interaction be- 
tween DM and baryons can also change the cosmic mi- 
crowave background (CMB) anisotropy and large scale 
structure (LSS) power spectrum. The actual constraints 
depend on the coherent or incoherent properties of the 
scattering. Taking the more stringent limits for a coher- 
ent scattering (o~ X He — 16c X p)' the u PP er constraints at 
95% C.L. are shown in the purple line of Fig. [1] 20]. 

Bound States: If SIMPs can form bound states with 
nucleons, then these states can be found in heavy water 
searches. These constraints are very severe, and would 
exclude such models. We must make sure that such 
bound states do not form. 

The simplest way to make this happen would be if the 
force between DM and nucleons is repulsive. This could 
happen if the force mediator was a vector, and the DM 
and SM particles all had the same sign of charge. In 
that case, however, the charge neutrality of the universe 
is violated, and while this is not obviously impossible, it 
is difficult to find models where a nonzero gauge charge 
can be generated. 

The alternative is to take the mediator to be a scalar 
and the nucleons and DM to have opposite signs of 
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FIG. 1: Current experimental limits on SIMPs (see text for 
details). The regions above the black and purple solid lines 
and the regions enclosed by the blue and green contours are 
excluded. 



the couplings. We are therefore led to a model with a 
fermionic DM particle of mass m x interacting with nu- 
cleons through the exchange of a mediating scalar <f> of 
mass m<p. The interaction Lagrangian is taken to be 



£mi = - 9 X 4>XX ~ 9n 4> NN , 



(1) 



where N — p,n. For positive m x and and g x gN < 0, 
the force between DM and nucleons is repulsive at tree 
level, and no bound states form (we assume that the loop- 
level generated attractive force is subdominant). 

One also needs to consider additional constraints on 
the DM self-interacting strength from halo shapes and 
merging galaxy clusters such as the Bullet Clusters. This 
constraint sets (22, 23] 



< 3 GeV" 



(2) 



which also ensures that the DM does not self-scatter suf- 
ficiently many times to form bound states, even though 
the force between them is attractive. 

Since the scalar mediates an attractive force between 
nucleons, they can in principle form bound states as well. 
This is more complicated to analyze because there is al- 
ready a force between nucleons due to pion exchange 
which cannot be calculated precisely. We will therefore 
require the new scalar force to be small compared to two- 
pion exchange; we conservatively take m$ ~ m K and 
9n Ss 0.3 g-nNN which would lead to small corrections to 
the nucleon potential [24| ■ 

For the DM mass equal or below the proton mass and 
from the constraint in Eq. ((2|), we find 



< 



<? 2 



1 GeV 



10" 27 cm 2 . 



(3) 



where to calculate a xx and a xv we have neglected the 
non-relativistic DM momenta which are below m,j>. This 
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constraint is model-dependent; for g^ ^ g x an d m x < 
0.2 GeV one has even a stronger constraint than the ex- 
isting ones in Fig. [T] from the merging galaxy cluster con- 
straints. 

Summary: As can be seen from Fig. [TJ these com- 
bined constraints exclude SIMPs above 2 GeV (cf. 15]). 
On the other hand, for a low mass m x < 1 GeV, the 
constraints become much weaker due to the low-energy 
threshold limitation of direct detection experiments. Af- 
ter imposing all constraints, we still have a large region 
of parameter space allowed for SIMPs below 1 GeV with 
a scattering cross section as large as 10 -25 cm 2 . 

Collider Searches for SIMPs. The high-energy 
colliders do not have the limitation of direct detection 
experiments on searching for light DM particles. We now 
explore the possibility of seeing SIMPs at colliders and 
discuss their spectacular collider signatures. 

The first point to note is that a strongly interacting 
DM particle, if produced at colliders, may not be invisi- 
ble to detectors. Indeed, if the DM has interaction cross 
sections similar to hadrons, it would deposit energy in 
calorimeters, in a manner similar to neutrons and other 
neutral hadrons. More precisely, the signal depends on 
the scattering length L x , which is inversely proportional 
to the scattering cross section, of DM inside a detector. 
If, as is usually the case for WIMPs, L x is significantly 
larger than the calorimeter size, the DM leaves the col- 
lider without any trace and behaves as missing energy. 
However, for SIMPs with a large scattering cross section, 
L x could be smaller than the calorimeter size. The DM 
will then deposit most or all of its energy in the calorime- 
ters, especially in the hadron calorimeter (HCAL) be- 
cause it contains heavier elements than the electromag- 
netic calorimeter (ECAL). Therefore, we anticipate that 
the SIMP behaves like a jet at colliders. 

For a fast moving DM, its inelastic scattering cross sec- 
tion with nucleons (usually much larger than the elastic 
scattering cross section) determines the scattering length 
L x . We can relate it to the scattering length of neutrons 
inside materials via 

_incla 



A typical nucleon-nucleon scattering cross section is 
about 40 mb = 4 x 10~ 26 cm 2 . The bounds on the DM- 
nucleon elastic scattering cross sections, as noted above, 
are about 10~ 25 cm 2 , and the bounds should be much 
weaker for DM-nucleon inelastic scattering. It is there- 
fore indeed possible to have DM scattering lengths of sim- 
ilar order to hadronic scattering lengths. In fact, since a 
typical hadron scatters about 10 times inside the HCAL, 
the DM may appear as a jet even if its cross section 
is somewhat smaller than hadronic cross sections. For 
10~ 26 cm 2 < a^ u , the SIMPs will appear as jets. 

For the rest of this paper, we will focus on the region 
where the scattering length of the DM particles is much 
smaller than the calorimeter size, so that all the DM 
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FIG. 2: The DM jet signal production cross section for g x g u = 
1 for the solid and red line. The measured dijet invariant mass 
spectrum is shown in the circle and black points 27]. The 
squared and blue points are dijet invariant mass spectrum 
after applying the no-track and less-EMF cut on both jets. 



energy is deposited inside the detectors [31|. For pair- 
produced DM particles, we therefore have a dijet signa- 
ture. Furthermore, since DM is neutral, the DM jets 
should have no tracks and less electromagnetic activities 
in the ECAL and hence are trackless jets. Exactly be- 
cause of this property, the DM jets can be distinguished 
from QCD jets. 

To determine the feasibility of detecting these events, 
we first need to know the number of such events. This de- 
pends on further details of the model, in particular the in- 
teraction of the mediator with quarks and gluons (rather 
than nucleons). Using the matrix element of quarks in- 
side nucleons, one can relate the mediator coupling to 
quarks g q to the couplings to nucleons gw- For example, 
one has g^ = 8.22g u for N = p [2||. In terms of the 
coupling to quarks, we have the production cross section 
of <j(uu —> XX) — flxflu/(127rs) (where we have neglected 
the light DM mass). Using the MSTW parton distribu- 
tion function [26j . we have the signal production cross 
section at the 7 TeV LHC shown in the solid (red) line 

of Fig. m 

The background, which comes from QCD dijet pro- 
duction, has been measured [27j and found to be in good 
agreement with the SM prediction. We show the ob- 
served data in the circle dots (blue) of Fig. [2] as an esti- 
mate of the background. 

We have two main handles to reduce the QCD back- 
grounds. First, we require that there are no tracks in 
the tracker system. The charged particle multiplicity has 
been measured by the CMS collaboration at 0.9, 2.36 and 
7 TeV [28| and approximately follows the Koba-Nielsen- 
Olesen (KNO) scaling P(n) = (n)"^""/^ [29]. The 
average charged-particle number is (n) « 20 [28[, and 
hence approximately 5% of QCD jets have no charged 
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particles in them. Requiring two trackless jets therefore 
reduces the background by a factor of 400. 

We also require a small electromagnetic fraction 
(EMF) for each jet. For QCD jets, the EMF (coming 
from 7T° — > 77 and radiations of charged hadrons) is typ- 
ically peaked at 50% [SO] while for DM jets the EMF is 
proportional to the scattering length of DM in the ECAL 
over that in the HCAL and is around 10%. We can hence 
impose a cut that the EMF should be less than 20%. This 
further reduces the background by a factor of about 100. 
We note that there could be a correlation between the 
no-track cut and the EMF cut and the actual reduction 
on the QCD backgrounds requires detailed detector sim- 
ulations. 

Applying those two cuts, the QCD backgrounds can be 
reduced by a factor of around 40000; the reduced back- 
grounds are shown in the squared points (blue) in Fig. [2] 
We see that these cuts reduce the background to well be- 
low the DM production rate. Hence such DM jets can be 
discovered even at the early LHC. 

Conclusions. We have argued that light DM parti- 



cles which have strong interactions with the SM particles 
are not excluded by current astrophysical constraints. 
These particles can even have interaction cross sections 
close to hadronic interactions. We have also shown that 
these DM particles can have unique signatures at hadron 
colliders. Unlike WIMPs, these DM particles do not man- 
ifest themselves as missing energy; instead there are jets 
associated with them. We have proposed a new search 
strategy for SIMPs by looking for a signal of two track- 
less jets. We have shown that suitable cuts can plausibly 
reduce the QCD dijet background so that the DM jets 
may be discovered even in the early data of the LHC. 
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